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Abstract. The compactness of weak solutions to the magnetohydrodynamic equa- 
tions for the viscous, compressible, heat conducting fluids is considered in both the 
three-dimensional space R"^ and the three-dimensional periodic domains. The viscosi- 
ties, the heat conductivity as well as the magnetic coefficient are allowed to depend on 
the density, and may vanish on the vacuum. This paper provides a different idea from 
|15l to show the compactness of solutions of viscous, compressible, heat conducting 
magnetohydrodynamic flows, derives a new entropy identity, and shows that the limit 
of a sequence of weak solutions is still a weak solution to the compressible magnetohy- 
drodynamic equations. 



1. Introduction 

Magnctohydrodynamics (MHD) is the theory of the macroscopic interaction of electri- 
cally conducting fluids with magnetic fields. It has a very broad range of applications. It 
is of importance in connection with many engineering problems, such as sustained plasma 
confinement for controlled thermonuclear fusion, liquid-metal cooling of nuclear reactors, 
and electromagnetic casting of metals. It also finds applications in geophysics and astron- 
omy, where one prominent example is the so-called dynamo problem, that is, the question 
of the origin of the Earth's magnetic field in its liquid metal core. 

Due to their practical relevance, MHD problems have long been the subject of intense 
cross-disciplinary research, but except for relatively simplified special cases, the rigorous 
mathematical and numerical analysis of such problems remains open. In the viscous in- 
compressible case, MHD flow is governed by the Navier-Stokes equations and the Maxwell 
equations of the magnetic field. For the mathematical analysis in incompressible MHD 
equations, see [H [121 ES] and the references therein. In the compressible case, the math- 
ematical analysis is much more complicated, due to the oscillation of the density, the 
concentration of the temperature, and the coupling interaction of hydrodynamics with the 
magnetic field. The full system of the three-dimensional magnetohydrodynamic equations 
in the Eulerian coordinates can be read as follows ( [T71 II 8j ): 

pt -I- div(pu) = 0, (1.1a) 

{pu)t + div (pu (g) u) + Vp = (V X H) X H + div^-, (1-lb) 

£t + div(u(f ' + p)) = div((u X H) X H + z/H X (V X H) + u* + kV6'), (1.1c) 

Ht - V X (u X H) = -V X (i/V X H), divH = 0, (l.ld) 

where * = 2^iD{u) -f Adivul with 3A + 2p > and ^(u) = i(Vu + (Vu)^) denotes the 
strain rate tensor; p denotes the density, u e M'^ the velocity, H e M"^ the magnetic field. 
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and 9 the temperature; £ is the total energy given by 

E = p[e + l|u|2^ + i|Hp and E' ^p{e.+ \\n\'^ , 

with e the internal energy, |/o|up the kinetic energy, and ^|Hp the magnetic energy. The 
equations of state p = p{p, 6), e = e{p, 9) relate the pressure p and the internal energy e 
to the density and the temperature of the flow; I is the 3x3 identity matrix, and (Vu)^ 
is the transpose of the matrix Vu. v{p^ 9) is the magnetic field coefficient, k = 9) is 
the heat conductivity. In general, equations ()I.la|l . (jl.Ibp . ()I.lc|) denote the conservations 
of mass, momentum, and energy, respectively. The equation (|l.ld[) is called the induction 
equation, and the electric field can be written in terms of the magnetic field H and the 
velocity u, 

E = i/VxH-uxH. 

Although the electric field E does not appear in the MHD system (|l.la|) - (jl.ldp . it is indeed 
induced according to the above relation by the moving conductive flow in the magnetic 
field. 

In this paper, we are interested in the compactness of weak solutions to the compressible 
MHD equations (jl.ip both in the three-dimensional space and in the three-dimensional 
periodic domains. As it is well-known, the motivation of considering the compactness 
of weak solutions is to show the existence of weak solutions and the stability of weak 
solutions of nonlinear problems. In the literature, there have been a lot of studies on 
MHD by physicists and mathematicians because of its physical importance, complexity, 
rich phenomena, and mathematical challenges; see [il [51 [71 [^1 [T H [1 ^ [H I [T? l [TC I [T5 1 and 
the references cited therein. For instance, the smooth global solution near the constant 
state in one-dimensional case is investigated in [16j . However, many fundamental problems 
for compressible MHD with large, discontinuous initial data are still open. 

Positive results on the existence of weak solutions with large, discontinuous data for 
compressible MHD equations have been obtained recently in [T31 [H] , specially in [T3] for 
full compressible MHD equations with temperature-dependent viscosities. More precisely, 
it was shown in [15| . under certain structural hypotheses imposed on the pressure p and 
the heat conductivity coefficient k, that the full compressible MHD system admits at 
least a global-in-time variational solution for large initial data. Those solutions satisfy 
the equations p. lap , (jl.lbp , (|l.ld[) in the sense of distributions while the thermal energy 
equation p.lcp is being replaced by two inequalities to be accordance with the second law 
of thermodynamics. This approach is in the spirit of the concept of weak solutions with 
a defect measure introduced by several authors in different contexts, see [6]. However, 
in order to obtain the estimates on the gradient of the velocity, the works in [Ml [15] do 
rely strongly on the assumption that the shear viscosity p, is bounded below by a positive 
constant. 

Our aim, in this paper, is to show the compactness of weak solutions to the full com- 
pressible MHD equations with viscosity coefficients vanishing on the vacuum both in the 
three-dimensional space 'E? and in the three-dimensional periodic domains. Although the 
periodic case does not correspond to a physical configuration, its mathematical treatment is 
technically easier, while it retains the main mathematical difficulties of the problem of the 
flow. More importantly, in our context, the viscosities p, A, the heat conductivity k, and 
the magnetic coefficient v can be allowed to depend on the density p and the temperature 9 
of the flow. We remark that the similar problems for the compressible Navier-Stokes equa- 
tions have been studied in [H [21 [3l [21] . Comparing with those works on the compressible 
Navier-Stokes equations, we will encounter extra difficulties in studying the compressible 
MHD equations. More precisely, besides the possible oscillation of the density and the 
concentration of the temperature, the appearance of the magnetic field and the coupling 
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effect between the hydrodynamic flow and the magnetic field should also been taken into 
consideration. We remark that a simultaneous independent work similar to this paper for 
the case of the two-dimensional periodic domain was done in [22] . 

The novelty of this paper is to provide a new method to deal with the vanishing viscosities 
for compressible MHD flows. The loss of positivity of the viscosity coefficients implies that 
there is no hope to obtain directly the uniform bound on the gradient of the velocity. 
It is well known that the main difficulty, in proving the compactness of weak solutions 
of compressible MHD equations, is to pass to the limit for the nonlinear terms. In our 
context, the new kind of entropy equality will provide the estimates on the gradient of the 
density, which makes the nonlinear terms much easier to be dealt with and also give rise 
to a new estimate of pu^ in a functional space better than L°° {[0,T]; L^{^1)). In other 
words, although the regularity on the velocity that we can get directly is much lower, the 
regularity on the density in our context is much higher. To achieve this aim, the entropy 
equation (|3.9[) and the thermal equation (|4.2[) need to be taken into consideration. But, 
unfortunately, the case with constant viscosity coefficients is excluded from our setting and 
the extension to the general case in which the viscosity coefficients depend on both the 
density and the temperature seems also out of the reach of our present work. 

We organize the rest of this paper as follows. In Section 2, we will give the hypotheses 
in detail, introduce the definition of weak solutions, and state our compactness result 
(Theorem I2.1[) . In Section 3, we will derive the a priori estimates and a new kind of the 
entropy identity. In Section 4, some auxiliary integrability lemmas are showed. Finally, we 
will finish the proof of Theorem 12.11 in Section 5 using Aubin-Lions Lemma. 



2. Assumptions and the Main Result 

To our best knowledge, the rigorous mathematical analysis for compressible flows is 
beyond the available mathematical framework. Hence, we need add some restrictions to 
viscosity coefficients fi, A, the heat conductivity n and the magnetic coefficient v. 

2.1. Assumptions. To begin with, we assume that fi{p) and A(p) are two C^(0,oo) func- 
tions satisfying 

A(p) = 2(p^'(p)-Ai(/9)). (2.1) 

As seen later on, this relation is fundamental to get higher regularity on the density. More 
precisely, with the help of this relation, we can show a new kind of entropy equality, which 
then gives the uniform bound on the gradient of the density. Next, due to our technical 
restrictions, we will need the following constraints: /i(0) = and there exist positive 
constants cq, ci, A and m>l,2/3</3<l such that 

co_ 

CI 

Observing that the assumption (|2.2p implies that fJ,'{p) > for p > 0. 
The heat conductivity coefficient k is assumed to satisfy 

K{p,9)^Ko{p,9){p+lW + l), (2.3) 

where a > 2, and kq is a C°(i?^) function satisfying for all positive p and 0, 

C2 < Ko(p, 6*) < — , 

for some positive constant C2- 

For the pressure, we assume that the equations of state are of ideal polytropic gas type 

p^pe + p,{p), e^c,9 + Peip), (2.4) 



for aU s < A, cqs^-i < /i'(s) < ^ and 3A(s) + 2p{s) > cqs^^, 
for aU s > A, cis™"! < n'{s) < and cis" < 3A(s) + 2p{s) < 



(2.2) 
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with Pe{p) = />e(e)/e'dC. Wc also require that Pe{p) satisfies 

f C3p-'-i<p^(p)<^p-'-\ ifO<p<^o, ^2 5) 

\ P'eip) < cap"-', iip>Ao, 

for some > 0, C3 > 0, C4 > 0, Z > - 1, and k < (m - i) ^^^z^. 

For the magnetic coefficient, we need the following assumption: 

1 
J^{p, 0)>C5- on {p > 0}, - > i^ip, 9) > C6, (2.6) 
P ce 

for some C5 > 0, ce > 0. 

2.2. Main Result. Before we state the compactness result, we need to specify the def- 
inition of weak solutions which we will address. It is necessary to require that the weak 
solutions should satisfy the natural energy estimates and from the viewpoint of physics, 
the conservation laws on mass, momentum and energy also should be satisfied at least 
in the sense of distributions. Based on those considerations, the definition of reasonable 
global-in-time weak solutions goes as follows. 

Definition 2.1. A vector (p, u, 6',H) is said to be a global-in-timc weak solution to the 
full compressible MHD system (jl.lap - (jl.ldp . if and only if for any positive number T, the 
following conditions are satisfied: 



p£ e L°°{[0,T];L\n)), p\uf G {[0,T]; LH^)), 
e L°^([0,T];L^{n)), {p"^^ + p"'^^)Wu e L^{[0,T];L^(n)), 



{l + ^)\/9^/^eL'{[0,T];L'm, {1 + e L\[0,T]; L^nj), 

for a > 2. Moreover, for large enough s > 0, we have 

p, pu, p£, UeCi[0,T];H~'{n)), 
and hence, the initial data are satisfied in the sense of distributions (denoted by 

v'{n X (o,r)). 

• The equations (|l.lap - (|l.ld[) are satisfied in the sense of distributions. 

Now our compactness result can be read as follows: 

Theorem 2.1. Let be either the three-dimensional periodic domains or the three-dimensional 
space K^. Assume that p, A, v, k are C^[0, 00) functions satisfying the hypotheses (|2.1[) - 
(j2.6p . Let {(p„, u„, 0„, H„)}^]^ be a sequence of weak solutions of (jl.lap - (|l.ldp satisfying 
the entropy equation p.9p and the thermal energy equation (j4.2p with the initial data 

Pn{x, 0) = Po,n{x), Unix, 0) = Uo,„(x), 

6'„ (x, 0) = 6'o,„ (x) , H„ (x, 0) = Ho,„ (x) , 
where po,n, Uo,„, 6*0, „, Ho,„ satisfy 

Po,n>0, po,n ^ Pq in L^{Q), 

P ^PoluoP znLi(f^), (2.7) 
Oo.n > 0, Ho,„ ^ Ho m L\n), 
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Assume also that the sequence of the magnetie field {H„}5^j is uniformly bounded in 
L°°{n X {0,T)). Then, up to a subsequence, {(/9„, u„, H„)}J^j^ converges to a global- 
in-time weak solution (p, u, 0,H) with initial data (po, uq, Ho). More precisely, 

{Pn}^=i converges strongly to p in C{[0,T]; LP{il')) for all 1 < p < 6m— 3; 
{u„}^=i converges weakly to u m L9i([0,T]; W^i'?3(f7')) for qi > 5/3, 93 > 15/8; 
' {(^n}^=i converges strongly to 9 in LP([0, T]; i'^(ri')) for all p < a and q < Sa with 
- = + TT^, for all r < 3/2; 

q p{a—l)r 3p{a — 1) ' ^ / ' 

{H„}^^i converges weakly to U m L^{[0,T]; H\n)) Ci[0,T]; Ll,^^^{n)), 
where D,' is any sufficiently smooth and compact subset offl. 



3. Energy Estimates and the Entropy Inequality 

In this section, wc dedicate to the well-known a priori estimates and a new kind of 
the entropy equality on weak solutions of the compressible MHD system (|l.lal) - (jl.ldp . To 
begin with, from the total energy equation (jl.lcp . the physical energy inequality holds 



Eit) :=^^p(^e + i|up^ it,x) + ^\U\\t,x)dx 



(3.1) 



As shown in [Ml [15], the energy estimate (|3.1[) alone is not sufhcient to build up a 
reasonable compactness theory of weak solutions to compressible MHD equations in the 
sense of distributions since we can not obtain any a priori estimate on the dissipation about 
the viscous stress and the gradient of the magnetic field. Comparing with a priori estimates 
for isentropic cases (see [H]), this is a major difference, because in the isentropic case, the 
viscous dissipation naturally provides a bound in spatial variables on the velocity u. 
To establish the compactness theory in our new framework, the following calculation is 
crucial: 



Lemma 3.1. 

ld_ 

2di 



and 
1 d 
2di 



[ (p|u|2 + |H|2) dx+ [ 2p{p)D{u) : D{u)dx 
Jn Jn 

+ / X{p)\divu\'^ dx + / i^|VxHp(i.T= / p{p,9)diYudx, 
Jn Jn Jn 

[ (/9|u + 2V(/j(p)p + |H|2) dx+ [ 2p{p)A{u) : A{u) dx + [ x dx 
Jn Jn Jn 

[ p{p, 6')divudx- - 2 /" Vp{p, 9) ■ Vip{p) dx + 2 [ (V x H) x H • dx, 
Jn Jn Jn P 



(3.2) 



(3.3) 



where A{u) = (Vu — Vu^)/2 denotes the skew symmetric part o/Vu, and(p'{x) — p'{x)/x 
for all X > 0. The notation A : B denotes the dot product between two n x n matrices A 
and B . 



Proof. The energy equality (|3.2|) is classical, and can be shown by multiplying the mo- 
mentum equation ()l.lbp by u, the mass conservation equation (jl.lap by |up/2, and the 
magnetic equation (jl.ld[) by H, then summing them together. Here we used the following 
identity: 

( (V X H) X H • udx = - / (V X (u X H)) • Hdx. 

Jn Jn 
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Now, we turn to show the equality p.3p . The idea is taken from [T], and the argument 
goes as fohows. From the mass conservation equation, we deduce that 

dtipip) + u • Vip{p) + Lp'{p)pdYvu = 0. 

This gives, differentiating this equation with respect to the space variable Xi, i = 1,2,3, 
noting X = (xi, a;2, xa), 

dtd^(p{p) + (u • W)d,(p{p) + id,u ■ V)(/9(p) + a,(v7'(p)pdivu) = 0. 

Let us multiply this equation by pdiip{p) and sum over i, by using the mass equation, then 
one can deduce 



(3.5) 



(3.6) 



[ p\Vip{p)\'^ dx + [ pVifiip) i^Vtpip) -.Vudx + [ y{ip'{p)pdivu) ■Vp{p)dx = 0. 
2 dt Jq 

(3.4) 

Multiplying the momentum equation by \7fi(p)/p, we get 

/ (a.u + u . Vu) . Vpip) dx + 2f pip)Din) : ( - YM^] dx 
Jn Jn VP J 

+ / Vp{p,0)-^p'{p)dx + 2 f V{{ii{p)~ p,'{p)p)diYn)-^-^dx 
Jn P Jn P 

= /(VxH)xH.^dx. 

Jn P 

Integrating by parts, equation (j3.4[) can be rewritten under the form 

/ p|V^(p)Pdx+ / Mp)VMp) -Vu-Vp^^ f Mp)Vu : WMp) 
2dt7n Jn P^ Jn P 

- I Mp)Vdivu ■ VMp) ^ / V(v.'(p)pdivu).VMp)rfx = 0. 
Jji P Jn 

Adding equation (|3.5[) to equation (|3.6p multiplied by 2, we get 

/(9.u + u.Vu).VMp)rfx-2 / Mp)Vdivu.VMp) ^^ 
Jn Jn P 

+ ~ f 2p\\/^{p)\'dx + 2 I X/{{p{p)-p'{p)p)dWu)-^^dx 
2 dt Jq Jo P 

+ 2 / V(^'(p)divu) • VAi(p)da;+ / S/p{p,e) ■ —p.' (p) dx 
Jn Jn P 

^ /(VxH)xH.^d.. 
Jn P 

By splitting the terms involving divu and by summing them, we get 

/ {dtu + u-\/u) ■Vp{p)dx + ~ I 2p|Vv?(/9)pdx+ / Vp{p^e) - ^p' {p) dx 
Jn ^dtJa P 

= /(VxH)xH.^<ix. 
Jn P 

But as for the first term in (13.81). we can calculate 



(3.7) 



(3.8) 



f {dtU + u-Vu)-Vp{p)dx ~ [ VL-Vp{p)dx— [ VL-Vdtp{p) dx+ f {u-\7u) -V p{p) dx. 
Jn Jo Jn Jn 
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By using now the mass equation and integrating by parts the last two terms, this gives 
/ (dtU + u-Vu) ■V^i{p)dx = — I u-Vn{p)dx~ f /i'(p)div(/9u)divuda; 

~ / /i(p)u • Vdivuda; — / fi{p)diUjdjUi dx. 
Jn Jn 

Integrating by parts in the third term, we get 

/ (9tU + u • Vu) • V//(/9) =-^ / u-\/p{p)dx~ ( {pp! {p) — p.{p))\A\vu\'^ dx 
Jn Jn Jn 

- / p{p)diUjdjUidx. 
Jn 

Adding the above identity with (|3.8p . we get the foUowing equaUty 

/ u • Vp{p)dx - f (pp'ip) - p{p))\divu\'^dx - f p{p)d,UjdjU,dx 
Jn Jn Jn 

+ ~ I 2p\y^{p)\^dx+ f yp{p^9)-^p'{p)dx 
2 dt p 

= /(VxH)xH.^dx. 
Jn P 

Adding this last equation multiplied by 2 to the energy estimate p.2p gives p.3p . □ 
Next, we introduce the concept of the entropy s(p, 9) which satisfies the entropy equation 

dt{ps) + div(p.u) + div ( ^\ = i(HV X Hp + * : Vu) - ^^^^'7^'' , (3.9) 



with s(p, 6) = \n9 — Inp, where is a positive constant denoting the specific heat at 
constant volume. The entropy equation is useful in compressible flows, because it provides 
naturally the estimates in the gradient of the temperature. More precisely, integrating 
(|3.9p over x (0, i), the following proposition is verified: 



(3.10) 



Proposition 3.1. Assume that poso G L^(ri). Then, for all t > 0, one has: 

j'j^ i(i^|V X Hp + * : Vu) + ^ii^l^^dxdi < j^Ps+ \poso\dx. 

Observing that 

ps < CvpO — pin p. 
the first term on the right-hand side of (j3.10|l can be estimated by 



psdx < / Ci,p9dx — I plnpdx. 
! Jn Jn 

Multiplying (|l.lap by 1 + In p, we get 

dt{p\np) + div(pulnp) + pdivu = 0. 

Thus, we have 

plnpdx— / po\npodx+ / / pdivudxds. 



8 



XIANPENG HU AND DEHUA WANG 



Therefore, the right-hand side of (|3.10p can be estimated by 



ps dx < 



< 



Cupd dx 
C upd dx 



\Po In Po\dx 
|/9o Inpol^a^ 



p|divu| dxds 
V3A + 2p ^/O 



dYvvL\yJ p9 dxdt, 

(3.11) 

and, then using the classical Young's inequality, the bound of p9 in L°°([0, T]; 

and the assumption (|2.2p that ensures that s ^ s/3A(s) + 2/x(s) belongs to L°°{R+), we 

conclude that the terms on the left-hand side of (|3.10p is bounded in L^{n x (0,T)). 

Hence, if po^o f-nd polnpQ belong to L^(f2), then the components of the following four 
quantities ^/3X + 2/i|divu|/v^, ^£)(u)/v^, V^^V x H/v^, {^/p + l)V6't and (^p -t- 
l)Vln6' are bounded in L'^{n x (0,r)). We note that the last two bounds involving the 
temperature gradient provide the following useful estimates: 

(Vp + l)V6l° e L^in X (0, T)), for all a such that < a < a/2. (3.12) 

In order to get enough a priori estimates from Lemma 13.11 and the initial condition 
p.7p . we have to control the following terms: 

Vm(p) 



pdivudx, Wp-Wip{p)dx, / (V x H) x H ■ 
Jn Jn Jn 

To this end, the following estimates are useful: 



■ dx. 



Lemma 3.2. Let fl be the three-dimensional periodic box or the whole space R'^. For all p 
satisfying p~^/^Vp(p) S L^{fl), one has 



\r-'^'X{p>2A}\\LHn)<c\ 

V3-1/2 



Vm(p) II 
-^\\l^(Q), 

vMp) I 



i il/-^/'X{p<A/2}llLa(0)<c||^|U.(a), 

for some positive constant c, where A is from p.2|) . and x is the characteristic function. 

Proof. Let us consider the function r/ = a^™"^/^, where ^ : (0,oo) — > (0, cx)) is a smooth 
increasing function such that ^(s) = s for s > 2A and ^(s) = for s < ^ and a is a positive 
constant. By hypothesis (|2.2p . we can choose c such that ri'{s) < cfj,'{s)/y/s for all s > 0. 
Using Sobolev's inequality, we have 

Vp(p) 



|f?(p)||Le(ji) < c|lV?7(p)|li2(j2) < c 



VP 



L2(0) 



The left-hand side of above inequality is bigger than j|p™ ^^^X{p>2A}\\L<i{n)- This imphes 

V/i(p) 



l|p'"-'/'X{p>2A}llLe(^2)<C 



VP 



To prove the second part, a similar approach can be applied. Indeed, choosing the 
function 77 = a£,^~^^^ such that |77'(s)| < cp'{s)/^/s for all s > 0, where ^ : (0,oo) (0,oo) 
is a smooth positive function such that ^(s) = for s > A and ^(s) = s for s < A/2. By 
Sobolev's inequality, we have 

Hp)\\LHn)<c\\\/vip)\\LHn)<c ^^^^^ 



L2(0) 

The left-hand side of above inequality is bigger than \\p^~^^^X{p<A/2}\\L<^{n)- This implies 



1/ ^^^X{p<A/2}llL6(n) < c 



VP 

,/3-l/2 

Vm(p) 



VP 



L2(n) 
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□ 



From Lcmma[33 we know p G L°°([0, T]; for any bounded subset n' of n. 

Lemma 3.3 (The control of J^pdivudx). 



pdivudx < 



dt 



I pPe{p) dx + eIIa/SA + 2fidh 



L2 



l2 



l2 

L6 



l2 



Vp{p) 



VP 



for all positive e. 

Proof. By the continuity equation in the renormahzed sense, we have 

/ pdivudx= / pe{p)divudx + / pOdivudx 
Jn Jn Jn 



d 
It 



pPeip) dx 



(3.13) 



pOdiyudx. 



For the second term on the right-hand side of (|3.13p . we can estimate 



p^divu dx 



< ||v/3A + 2/idivu|| 



L2 



X {\\pOX{p<a}/V3^ + M\l- + \\p0X{p>a}/V^^ + Ml-) 

< c|| ■\/3A + 2/idivu||^2 

X {\\py'0U4p^'^-'^^y'\y^A}\\L^ + A-^--/^pX^^^^^^ 

< c|| ■\/3A + 2/_tdivuj|i2 

X {\\p9\\l{'\\0\\l^U('-^^^^^^^ 
Thus, in view of Lemma 13.21 and Young's incquahty, we have 

< c|| ■\/3A + 2/xdivu||^2 



pOdivu dx 



X (\\p0\\l^ + \Ml<^ + 
< e|| ■\/3A + 2^divuj||2 



Vm(p) 



LI 



2 



2 

L3 



VP 



Vm(p) 



V~P 



□ 



Bv (|3l2l) and Sobolev's inequality, 9 e L^{[0,T]; L^{n)) n L^{[0,T]; L^{n)), since a > 2. 
Also, in view of the total energy conservation inequality p.ip and our assumptions, we 
know that pPeip) e L°^{[0,T];L^{n)) and p9 G L°° ([0,T]; L^{n)). Thus, it is possible, by 
taking e small enough, to get some a priori estimates from p.2p and ()3.3p via Gronwall's 
inequality. 



Lemma 3.4 (The control of j^Vp ■Vif{p) dx). 

- [ \Ip-\/ip{p)dx <- [ ip' {p)e\S/ p\^ dx + [ [ceK{p,9) 
Jn Jn Jn \ 

1 + 1-/3 „ 

|Vp 2 \^X{p<Ai}dx, 



dx 
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with Ai = min{^, Aq}. 

Proof. By our assumption (|2.4[) . we have 

Hence, we have 



/ Vp ■Vip{p)dx = - / if' {p)e\V p\'^ dx - / Lp'{p)pS/e - Vpdx 
Jn Jn Jn 

' I p'ApW[p)\yp?dx 

< - [ ip'{p)0\\/p\^ dx - ( ip'{p)p\/e -S/pdx 



(3.14) 



c / |Vp- 
Jn 



'X{p<Ai} 



dx. 



because f'{p) > and p'^ip) > 0. 

As for the second term on the right-hand side of (j3.14p . wc have, by our assumption 



ip'{p)p\/0 ■ \/pdx 



< [ \ip'{p)p\ie ■\ip\dx 

Jn 



< ^ (^CeK{p,( 



\ve\^ , pe^ \vp{p)\^^ 



Kip, 9) p 



dx 



da:. 



Thus, we have 



/ Wp-\/ip{p)dx / ip' {p)9\\7 pI"^ dx + / [CeK{p. 
Jn Jn Jn \ 



dx 



1+1-13 ,o 

|Vp ^ \ X{p<Ai} 



dx. 



□ 



Noting that Proposition 13. II imphcs n{p,6)^-^^- G L^{Q x (0,r)). Therefore, it is also 
possible, by incorporating the estimate into p.2p and p.3p . to get some a priori estimates 
via Gronwall's inequality. 

Lemma 3.5 (The control of ^(V x H) x H • dx). 



( (V X H) X H 

Jn 



V/i(p) 



dx 



< c 



|VxH|2t.(p,^,H) , |V/i(p)|^ 



dx. 



Proof. Indeed, we can estimate, by our assumption (j2.6p and the uniform bound of H„ in 
L°°{n X (0,T)), 



/ (V X H) X H • dx < f |(V X H) X H 

Jn P Jn 



V/^(p) , 



dx 



< c 

< c 



|VxH|2|H|2 , |Vm(p)|^ 



dx 



P P 

IVxhiMp,^) I |vMp)P 

9 ^ p 



dx. 



□ 
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The entropy inequality ((3l0)) implies that l'^>^H|Jt/(p,e) ^^i^^^^^ ^.^(l^ x (0,r)), pro- 
vided pqSq and polnpo belong to L^{fl). Therefore, from Lemma 13.1113. 51 we can deduce 
the following a priori estimates via Gronwall's inequality: 

\/puj|i,= ([o,T];L2(i2)) < C, ||p"^/^V/l(p)||ioo([o,T];L2(J2)) < C, 
|(p™/2+//2)Vu|U2(ox(0.T)) <C, r^'y Vu|U2(ox(0,T)) <C, 
(0)l/2(^(/3-l)/2+^(™-l)/2)^||^^^^^^^_^^^ < ||pPe(p)|U=e([o.T];Li(n)) < C, 

y^^V(^(p)||i,2(nx(o,T)) < c, !|p6'||L~([o,T];Li(a)) < c, 



H||L=o([o,T];L2(n)) < C, ||\A'V X H||i2(nx(0,T)) < c, 

2(nx(o.T)) < c, 

(3.15) 



= ([0,T];L2(n)) S C, IIV^^V X ±l||i2(nx(0,T)) 
a + pV6'"||i2(j^x(0,T)) < C, l|Vp ^X{p<Ai}||l2(Ox(0,T)) < c, 



for all a £ [0,a/2]. 

4. Some Integrability Lemmas 



As mentioned in [151 120j , the lack of a priori estimates on approximation solutions to 
the compressible flow is the main difficulty to prove the existence and the compactness of 
global-in-time weak solutions. Indeed, the basic and natural a priori estimates are not 
sufficient, since the energy equation does not hold so far even in the distribution theory 
framework. For more details, we refer the readers to [15] . 

This difhculty has been circumvented in [ini [TS] by restricting the generality of the 
equations of state (|2.3p and (|2.4|) . and defining variational solutions for which the energy 
equation ([l.lcp becomes two inequalities in the sense of distributions. However, this ap- 
proach requires significant restrictions on the equations of state, in particular the ideal gas 
case is excluded. 

This section is devoted to the local integrability analysis of the various energy fluxes such 
as pu|up, pue, up, kWO. One of the crucial steps is the additional integrability obtained 
on p. 

4.1. Integrability of the velocity. Let us begin with some bounds on the velocity with 
density dependent weights. 

Lemma 4.1. Let O be either the whole space M.^ or the three-dimensional periodic box, 
and T > 0. Let u be a vector field over n x (0,T) such that u e X«i ([0, T]; Lg^(r2)), 
ypu e L~([0,T];L2(f2))^ and p€ L°°{[0,T];LP^^{n)) such that 

, X 1 2(71 

gi e 1, 2 , and - + , ^, < 1. 

p q2{qi ~ 1) 

Then, there exists S > 3 such that p^^^u <E L^{fl' x (0, T)) for all bounded subsets Q' in fl. 

Proof. For the proof we refer the reader to Lemma 6.1 in [1]. □ 

In order to apply Lemma 14.11 to improve the integrability of p^^^u, we need first to 
deduce the integrability of the velocity u. Indeed, following the computation in [I], one 
may write Vu = p^^^"^ p^/'^Vn, and then deduce that 

ll'^u||^gi([0^T]:L<J3(O')) 

< C{n') {l + !|p"''/^X{p<Ai}IIl2j([0,T];L«(O'))) l|P^^^Vu||i2(f2x(0,T)) (4.1) 

< C{n') (l + \\Vp-^' + '-^^'\{p<A,}\\L-i[W{n'))) IIP''/'Vu|U2(Ox(0,T)), 



with 



'^'-^ .1=2^1- and ' ' ' ' 



/3 V ^ + 1/ 93 6j 2 
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In Lemma l4.ll letting (72 = ^ taking p = 6m — 3, we deduce that p^^'^u G 

L''((0,T) X r^') for some (5 > 3. 

As a byproduct of previous analysis, we also can derive useful bounds on energy flux. 
More precisely, one has for all bounded subset fl' of fi, 

u|U = ([0,Tl;L'-(n')) < C'IIp" llL-([0,rl;Li(O'))"^" \\L^:ii[lO,T]-L<^^i{n'))\\^\\L-i{[0,T]-L''2(^Q')), 

with 

._ 1 + 1-13 1_ 51 + 3 1 _ 17Z + 15 - 12/3 

^ 21 ' s ~ 6(« + l-/3)' 7 ^ 18(1 + 1- (3) ■ 

Noting that the hypothesis (|2.5p implies that s > 1 and r > 1. 

In order to bound the energy fluxes, it remains to control p'^u in L^{{0, T) x 17') for some 

5>1. Since u is bounded in L«i ([0, T]; L«^"^(r2)) and p is bounded in L°°([0, T]; 

the hypothesis (|2.5p implies that p^u is bounded in L''((0,T) x f^') for some 6 > 1. 

4.2. Integrability of the heat flux. In this subsection, we will need the following inte- 
grability on the temperature: 

Lemma 4.2. Let Q, be either the whole space or a three-dimensional periodic box and 
let T > 0. Let be a function over n x (0,T) such that {^+1)V0°-^'^ and (^+ 1)V In 61 
belong to L^{n x (0,T)) with a > 2, pe e L°°{[0,T];L\n)), and p-^/'^Vp.{p) G L^{n). 
Then, 6^"^ belongs to L^{[0,T]; L^{Q)) for allO<c<l. 

Proof. For the proof, we refer the reader to Lemma 7.3 in [T]. □ 

At this stage, in order to improve the integrability of the heat flux, we need to use the 
following thermal energy equation (cf. equation (1.15) in [15]). 

dt{p0) + div(p6lu) - div(K(p, 0)V0) = i/|V x Hp + : Vu - 6lpdivu. (4.2) 

As a matter of fact, we have 

Lemma 4.3. For any nondecreasing concave function from i?+ to R, one has 

[ f{0){2fiD{u) : D{u) + \\divu\^ + i^\S/ xli\^)dx- [ K{p,0)f"{0)\\/0\^ dx 
Jn Jn 

<^f pf{0)dx+ f p0f{0)\divu\dx. 
"I Jn Jn 

Proof Multiplying (g^l) by f'{0), one has 

[ f{0){2^iD{u):D{u) + \\divuf + i'\Vxlif)dx- [ n{p,0)f"i0)\V0f dx 
Jn Jn 

f'{0){0pdivu + dt{p0) + div(p6'u)) dx 



f'{0){0pdivu + pdt0 + puV0) dx 

f'{0)0pdiYU + pdtf{0) + pnVf{0) dx 

f{0)0pdivu + dt{pf{0)) + d\Y{puf{0)) dx 

< f f{0)0p\divu\dx + ^ f pf{0)dx, 
Jn at Jn 

here, we used twice the mass conservation equation (jl.lap . □ 
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Now, wc consider f'{9) ~ for some < c < 1 in Lemma l473l then we have 

^^{2^iD{n) : D{n) + Ajdivup + x Hp) dx + ^)^|V^?P rfa; 



< 



dt I P 



dx 



Idivul da;. 



Keeping the hypothesis (j2.3p in mind, we have 



(l + p)|V(l + 0)('^-^+i)/2|2rfa;< 



(1-c) 



dx+ / pe^-''\d:iY\i\dx. (4.3) 



For the second term on the right-hand side of (|4.3p . one has 
p^^-'^divu dx 



< ||v/3A + 2^divu|li,2 

X {\\pO'-'x{p<A}/V3X + 2ii\\L2 + ||p0'-^X{p>A}/v/3A + 2^|UO 

< c\\y/3X + 2^divu||i2 

< c\\^/3X + 2^divu||i2 

Thus, in view of Lemma 13.21 and Young's inequahty, we have 



p9^ '^divudx 



< c||V3A + 2/.divu|U2 + II^'^^IIl^ + 



< e||v/3A + 2^divu| 



h+c. 



V//(p) 


' ) 


VP 





Because pd and p belong to L°°{[0,T]; L^{fl)), we deduce that pO^"^ ~ p'^{p6Y~'^ belongs 
to L°°([0, T]; L^(r2)). Hence, integrating over t in the both sides of (|4.3p . and combining 
the estimates in Proposition 13. 11 one has, if < c < 1, 



(1 + p)|V(l + 0)(«-^+l)/2|2 dxdt < C, 



(4.4) 



"'SI 



for some positive constant C. In particular, by Sobolev's inequality, 9°^ belongs to 
i2([0,T];L3(i7)). 

Lemma 4.4. Let O be either the whole space M.^ or the three-dimensional periodic box and 
let T > 0. Assume that 9 satisfies the conditions in Lemma \4-2\ and yT~+~pV0'°~^+^-'/^ G 
L'^{[0,T];L'^{n)), 9"-"+^ e L^{[Q,T]-L^{n)), p € L°°([0, T]; L6"i-3(^)) y/jg„^ one has, 
for some p > 1, 

K{p,9)V9(^L^{[Q,T]-L^{n)). (4.5) 



Proof. For the proof, we refer the reader to Lemma 6.2 in [T]. 



□ 
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5. Compactness of Weak Solutions 

With the a priori estimates and integrabihty lemmas obtained in the previous sections, 
we now study the compactness of sequences of weak solutions {(p„, u„, H„)}^]^ and 
pass to the limit in nonlinear terms. 

To begin with, we will state the Aubin-Lions compactness lemma fseejl9j. Ch. IV, and 
[21] for more recent references) which we will use later. A simple statement goes as follows: 

Lemma 5.1 (Aubin-Lions Lemma). Let T > 0, p £ (l,oo) and let {/„}^i be a 
bounded sequence of functions in LP([0,T]; H) where H is a Banach space. //{/nj^JLi 
also bounded in LP{[0,T];V), where V is compactly imbedded in H and {dfn/dt}^^i is 
bounded in LP{[Q,T];Y) uniformly where H C Y, then is relatively compact in 

LP{[0,T];H). 

5.1. Compactness of the density. From the uniform estimates derived in Lemma [3.21 
we deduce that the sequence pn is uniformly bounded in i°°([0, T]; i^"'^'^(f2')) for all 
bounded subset Q' of Q. Up to a subsequence, one may assume that p„ converges weakly 
to some p in L^([0, T]; Lf^^{n)). In fact, we have 

Lemma 5.2. 

dtipl') is bounded in L\[0,T]; L^^^{n')), 
Vp"^ is bounded m L°°{[0,T];L^/^{Vl')); 

as a consequence, up to a subsequence, pn converges almost everywhere and strongly to 
some element p in C{[0,T]; U'{VL')) for all 1 < p < 6m — 3. Moreover, the continuity 
equation (|l.lap hold in the sense of distributions. 

Proof. Let us consider the renormalized mass equation satisfied by h{p) = p™, 
df {h{pn)) + div(/i(p„)u„) + (m - l)p™divu„ = 0. 

The uniform bounds of in L°° {[0,T]; L^{n)) and ofpn^^^h{pn) in L°°{[Q,T]; L^{fl')) 

imply that /i(/9„)u„ is bounded in L°°{[0,T]; L^^'^{n')). On the other hand, p™/Mivu„ is 
bounded in L^{nx (0, T)), and the sequence p™/^ is bounded in L°°{[0, T]; L6(2m-i)/m(j7'))^ 
hence, p™divu„ is bounded in L'^{[0,T]; L^/^{^')). Thus, 9t(/i(p„)) is uniformly bounded 
inL2([0,r];L3/2(f7')). 

The spatial regularity of p™ can be estimated as follows. Since V(p™) = ^^^^fp^, thus, 

V(p™) is bounded in L°°{[0,T];L^/'^{n')). The estimates deduced above, and thanks to 
Aubin-Lions Lemma, give the strong convergence of in L-^([0, T]; L^(ri')). We denote 
the limit of by p™. Then, since the function h{s) is strictly increasing, we conclude that 
Pn converges strongly to p := (p'^y/"^ in L^([0, T]; L^(ri')), and hence, by interpolation, 
in L9([0,T];LP(f]')) for all 1 < p < 6m - 3 with i = i + (g^E^y^. And using the 
continuity equation l|l.lap again, we know actually p„ converges strongly to p (p™)^/™ 
in C{[0,T];LP{n')) for all 1 < p < 6m - S. 

Finally, we already know that {un}^^i is uniformly bounded in L'i^{[0,T];W^''^^{Q')). 
Thus, up to a subsequence, u„ converges weakly to some element u in L"^^ {[0,T]; VF^'*' (^^'))- 
By Sobolcv's compact imbedding theorem, we also have that u„ is uniformly bounded in 
L^^^([0,T]; L^(Q')) since qi > 5/3 and 53 > 15/8. As a consequence, we pass to the limit 
in the mass conservation equation: 

atp + div(pu) 0, inV'in X {0,T)). (5.1) 



□ 
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In the spirit of Lemma 15.21 we can pass to the hmit in the sense of distributions for the 
terms pnPeiPn), Pe{pn), PnPeiPn)un, Pe(pn)u„ since u„ converges weakly to some element 



5.2. Compactness of the momentum. In this subsection, we show the compactness 
of the momentum /9nU„ by Aubin-Lions Lemma. We already know from the previous 
subsection that pnU„ converges weakly to pu in L^([0, T]; L'^/^(r2')), due to the facts 



%:un G L°°([0,T];L2((7)) and p„ e L°°{[0,T]; L^(n')). To this end, wc need to es- 



tablish the uniform bound of dt{pn'^n) in some suitable functional space. Indeed, we can 
show that the sequence 9t(p„u„) is uniformly bounded in LP{[Q, T]; 77~^(ri)) for somep > 1 
and s large enough. 

From the momentum conservation equation (jl.lbp . we have 

9t(p„u„) = -div (p„u„ (E) u„) - Vp„ + (V x H„) x H„ + div'i'n- 

For the first term on the right-hand side, pnU„ (g) u„ is bounded in L^^^{[0,T]\ L^^'^{Q')) 
uniformly as a product of pl/^Un, bounded in L'^([0, T]; L'^(rj')) and p„, bounded in 
i°°([0,T]; L^™~^(ri')). For the second term, p„ ~ pipn,dn) is bounded uniformly in 
L°°([0,r];Li(f7)), since p„6l„ is uniformly bounded in L°° {[0 , T]; L'^ (fl)) and Peip„) is 
bounded in L°°{[0,T]; L'^{n)). For the third term, by ^J^, (V x H„) x H„ is bounded 
in L^{[0,T]; L^/'^{Q)). As for the fourth term div*„, from the fact ^/p{pn)D{un) and 
^|A(p„)|divu„ are bounded in L^(rix (0, T)), and that ^ pipn) and \/\X{pn)\ are uniformly 
bounded in L°° {[0,T]; L^{n')), we deduce that div^„ is bounded in L^{[0,T]; L^^^{n')). 
Therefore, the sequence 9t(p„u„) is uniformly bounded in L^^^ {[0,T];W~^'^ {Q')) . 
Next, we have 

^i^Pn'^nj^ — Pn^i^nj ~^ ^nj^iPn 



Using the hypothesis ()2.2p and the estimate p„ E L°°([0, T]; i^™ ^(Jl')), one can deduce 
that „,^P^ e L°°{[0,T];L^{n')) and ^^r^ G L°°{nx{0,T)). Hence, from the estimates 

in (j3.15|l . we deduce that pnU„ e L'^{[0,T];W-^'^{n')). Thus, by Aubin-Lions Lemma, we 
deduce that p„u„ converges strongly to pu in L^/'^{[0,T]; LP{D,')) for all 1 < p < 3/2. 

As a conclusion, the product p„|u„p, converges strongly to p|up in L^{Q x (0, T)), since 
PnUn converges weakly to pu in L°°([0, T]; L^/^{n')), strongly to pu in L^^^{[0, T];LP{n')) 
for all 1 < p < 3/2 and u„ is bounded uniformly in L^/^{[Q,T]; L^{n')). Using the fact 
pl/^Un = pl/^UnX{pr,<6} + Pn^'^7iPn^^'^X{pr,>e}, Pn^^u„ is the sum of a Uniformly small 
term in L^{fl' x (0,T)) and another term converging to p^^^ux{p:^e} and then we deduce 
that pl/^Un converges strongly in L^{[0,T]; L^{fl')) to p^^^u. Finally, by the interpolation 
and the uniform bound of pl/^Un in L^{n' x (0, T)) for some (5 > 3, we conclude that 
Pn^'^Un converges strongly in L'^{Q' x (0, T)) to p^/^u. 

5.3. Compactness of the temperature. In this subsection, we want to derive compact- 
ness results for the energy and the temperature 0„. The first step is to derive uniform 
bounds in LP{[0,T]; H^''{il)) for some p > 1 and s large enough for the sequence dt{£n)- 
Indeed, we can rewrite the energy conservation equation (jl.lcp : 

dt{£n) = -div(u„(£,'j+p„)) + div((u„ x H„) x H„ + i/„H„ x (V x H„) + u„*„ + k„V6'„). 

For the first term of the right-hand side, we already know that p„u„|u„p is uniformly 
bounded in L''{[0,T]; L'^{n')) for some g > 1. Also, we already get the uniform boimds in 



u in L«i([0,r];VFi'«3(O')). 




(5.2) 
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L''{[0,T]; L'^{^1')) with some q > 1 for the terms p„'u„, p^u„ and K{pn,dn)Vdn in Section 
4. Next, the uniform bound of p„u„ in L°° {[0,T]; L^^^{fl')) and the uniform bound of 6'„ 
in L'^{[0,T];L^{n)) impUes that p„u„6'„ is bounded in ^'([O, T]; for some g > 1. 

Hence, p„u„e„ and p„u„ are bounded in L''{[0,T]; L'^{Q')) for some q> 1. 
For the viscous flux Un^I^n, we note the facts that 



a/ p.{pn)D{un) and ^|A(p„)|divu„ are bounded in L^{n x (0,T)), 



and 



^Mp^p-^^^ and v^AG^p;;i/3 bounded in {[0,T]; L'''^^"'-^^^^^"'-^\n')) 

due to the hypothesis (|2.2p and Lemma [3T2l hence in L°°([0, T]; and plJ^Un is 

bounded in L3(rjx (0, T)). Thus, the viscous fluxes are bounded in L^/^{[Q, T]; ^^^/"(fJ')). 

As for the terms related to the magnetic field, we have the fact H„ is bounded in 
L°°{[0,T];L^{n)) n L^{[0,T],H^{n)), hence, by interpolation, in LP([0, T]; where 
p> 5, q < 30/11, and 

11 2 

g ^ 2 ^ 3p' 

Thus, (u„ X H„) X H„ belongs to L'^([0,T]; L'^{^1')) for some <? > 1, since u„ is uniformly 
bounded in L^/^{[0,T]; L^{n')). By the bound of the magnetic field coefficient i'{p,9), we 
know i/„H„ X (V X H„) belongs to L2p/(2+p)([o, T]; L^i/^^+i\n')) for p > 5, 2 < g < 30/11, 
hence, in L«([0, T]; L'?(17')) for some q>l. 

With the bound of dt{£n) in mind, we can show the strong convergence of the term Pn^n- 
To this end, we will follow the argument in [T]. First, we note that the strong convergence 
of y/p^Un to y/pu in L"^ {[0,T]; L'^{il')) and the strong convergence of PnPeXPn) to pPe{p) 
in L''([0,T];Li(f]')) for aU r G (l,oo). Let us introduce 

^ = {/e4c(f^)ll|V/||L^(o) = 1} 

and a sequence 7^ of regularizing kernels given for instance by convolution operators such 
that the following basic properties hold: 

sup||/-r,./|i < ^, 

and for all compact subset Vl' C il, there exists Ck.n' such that for all / G /C, 



\TkJ\\ 



k,n' 



and Tkf e for all s > 0. 



Then, we can deduce that for any compact subset f2' C f2, one has 

C . 

/a'x(0,T) 



-(l|Pri^'ri|lL3(o/x(0,T)) + II P^ll L2(n/ x (0,T)) ) 



I da; 



O'x(0,T) 



2ll(|H„P 



2ll(P"|u»l' 



?i||l1(J2'x(0,T)) 

' P|un'7fc6'„|lii(f2/x(o,T)) 
(p„Pe(Pn) - P-Pe(/O))^fe6'n||Li(J2'x(0,T)) 

(£„ - £)Tk6n dxt 

O'x(0,T) 



(5.3) 
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Let US observe that the first term of the above right-hand side is bounded by 
C 

--j:{\\Pn\\L=°([0,T\-L^(n'))\\(^n\\v'-{[0,T\-L^{n')) + IIpII L-'([0.T]:L3(O')) II ^11 L2([0,T];L6(0')) )■ 

Therefore, given e > 0, there exists an integer fcg such that the preceding term is less 
than e/4 uniformly in n. Dealing with the second term is an easy task since On converges 
weakly to 9 in i^([0, T]; L®(ri')), so that for n large enough, the second term is estimated 
by e/4. But \\TkoOn\\L'^(n') is uniformly bounded in n, whereas H^, p„|u„p and PnPe{pn) 
converges strongly respectively to H^, p|up and pPe{p) in x (0,T)) respectively, so 

that the sum of the third and the fourth term is estimated by e/4 for large enough n. For 
the last term with k = ko, the uniform bound of dt{£n) in L^{[0, T]; H~^{U)) where p > 1 
implies that up to a subsequence, converges strongly to £ in C([0, T]; if ""(fi)), so that, 
for large enough n, the right-hand side of above inequality is less than c. It follows that 
^JJhiOn converges strongly in x (0,T)). 

On the other hand, we know the strong convergence of p„ to p^^/'^ in ^'^([0, T]; U'{Vt')) 
for p < 6 due to the Lebesgue's dominated convergence theorem, the estimate p.lSp 
and the strong convergence of the density. And then, we deduce that 6'„ converges 
to 9 in Li([0,r];L[„^(f7')) for aU r < 3/2. Recalling the uniform bound of On'^ in 
L2([0,T];L6(0')), we deduce that 6'„ converges strongly to 9 in LP([0, T]; L9(ri')) for all 
p < a, and q < 3a with 

1 a — p P ^ 1 

q p(a — l)r 3p(a — 1) ' 

for aU r < 3/2. 

5.4. Compactness of the magnetic field. The aim of this subsection is to show the 
compactness of the magnetic field. From Lemma 13. 1[ and the hypothesis (|2.6p , we deduce 
that H„ e L'^{[0,T];H^{n))r\L°°{[0,T];L'^{n)). Thus, we can assume that H„ converges 
weakly to some element H with divH = in L ^{[0,T]; H^{n)) n L°°{[Q,T];L'^{n)). 

On the other hand, from the equation (jl.ldp . we know that 

5tH„ = V X (u„ X H„) - V X (z/(p„, 6i„)V x H„). (5.4) 

For the first term on the right-hand side of (|5.4p . we deduce that 

u„xH„eL5/3([0,r];Lio/^(r!')) 

because of u„ e L^/^{[Q,T]- L'^in')) and H„ e L°^{[Q,T]: L'^{VL)). For the second term 
on the right-hand side of (|5.4p . we have i^(p„,6'„)V x H„ G i^(r2 x (0,T)) due to the 
hypothesis (EH). Hence, 9tH„ is bounded in L^/^{[Q,T];W-^^^^/'' {VL')). By Aubin-Lions 
Lemma, we deduce that H„ converges strongly to H in L^/'^{[{),T]]U\Vl')) for any 5 < 
p < 6. Furthermore, due to the uniform bound of H„ in L°°([0, T]; L^(ri)), by using the 
interpolation, one obtain that H„ converges strongly to H in L^'([0, T]; iy'(ri')) for some 
p > 5 and some q > 5/2. 

5.5. Proof of Theorem 12.11 To finish the proof of Theorem 12.11 we need to check that 
the limit functions p, u, 9, H are indeed the weak solutions, as defined in the introduction. 
We will complete this proof by several steps. 

Step 1: Convergence of the mass conservation equation. Let us start with the mass 
conservation equation (|l.la|) . since p„ converges strongly to p in C([0, T]; L^i^Tl')) for all 1 < 
p < 6m — 3, and u„ converges weakly to u in L'^^{[Q, T]; W'^'''^{fl)), we deduce that, by the 
Sobolev's compact imbedding theorem, p„u„ converges strongly to pu in L'"([0, T]; L^(ri')) 
for some r > 1. In particular, the mass conservation equation (jl.lap is satisfied in the 
sense of distributions. 

Step 2: Convergence of the momentum conservation equation. For the momen- 
tum conservation equation (jl.lbp . the strong convergence of /9nU„ and PnU„ ® u„ in 
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L^{[Q,T]; L^{n')) can ensure the passing to limit in the sense of distribution for the two 
corresponding term in the momentum conservation equation (|l.lb[) . On the other hand, 
since H„ converges weakly* to H in L°°{[0,T]; L'^{n))r]L^{[0,T]; H^{n)), this implies that 
the nonlinear term (V x H„) x H„ converges to (V x H) x H in the sense of distributions. 
As a product of p„ and 6'„, which respectively converge strongly in C([0, T]; L'^{n')) and in 
Z/^(ri x (0, T)), the term V{pnSn) converges to the limit V{p9) in the sense of distributions. 
The term Pe{p) is already done in view of the hypothesis (j2.5p and the strong convergence 
of p in C([0, T]; U'{Q!)) for all 1 < p < 6m — 3. Thus, we are left to show the convergence 
of the viscous flux. In fact, 

^l{pn)D{vin) = D{p{p^)n.,,) - i fvP^u„ ® YAeA + Y^BA ^ ^u„) . (5.5) 

Since '^^j^ converges strongly to in L°°([0, T]; L^(n')) and y/p^u„ converges strongly 

to y^u in i^([0, T]; L^(ri')), the fist term on the right-hand side of (|5.5[) converges to the 
corresponding term in the sense of distributions. The convergence of the second term on 
the right-hand side of (|5.5p in the sense of distributions can be shown by using the weak 
convergence of Pn^^^V ^.{pn) to p^^^'^V p{p) in L'^{[0, T]; L^{fl')) and the strong convergence 
y/p^Un in L^([0, T]; L^(r2')). For the bulk viscous term A(p„)divu„, by the assumption 
(|2.ip . it may be written in the renormalized sense: 

A(/9„)divu„ = -2{dtp{pn) + div(^(p„)u„)), 

which can be shown directly by the convergence of p„ and u„, and hence, the convergence 
in the sense of distributions for the momentum conservation equation is done. 

Step 3: Convergence of the energy conservation equation. The main difficulties in this 
step lie in the passage to the limit for the energy flux u{8' the heat flux kV0, the 

viscous term u^", and the nonlinear terms (u x H) x H, z^H x (V x H), because we already 
showed that £n converges strongly to £ in C{[Q,T];H^'^{^)) for some s > 0. 

For the energy flux p„u„0„, since y/p^u„ and ^/]hiQn converge strongly in L'^(p! x (0, T)) 
to y^u and respectively, /9„u„0„ converges strongly in L^i^l! x (0,r)) to pwd. For 
the energy flux p„u„|u„p, the strong convergence of pn^^^ in C([0, T]; i''(ri')) for all 
p < 6 implies that pn converges strongly to 

in C([0,T];L3(J]')). Hence, the term 
Pn^^^v^p^u„ converges strongly to p'^l^^w in -^■^([0, T]; L^/^(ri')), because of the strong 

convergence of in x (0, T)). And Lemma [4.11 implies that p^'^Un is uniformly 

bounded in L\Vi! x (0,T)) for some (5 > 3. This fact, combining with the interpolation 
inequality, gives the strong convergence of pu|up in L^{^ x (0,T)). The analysis at the 
end of Section 4 tells the strong convergence of p„u„Pe(/On) and u„p„ in L^iVl! x (0,T)) 
to puPe{p) and up, respectively. Thus, the energy flux u„(f4 + Pn) converges strongly to 
u{£' +p) in L^n' X (0,T)). 

The strong convergence of 9 in LP{[0, T]; L'(il')) for p < a and q < 3a implies that 9f/^ 
converges strongly in L^([0, T]; L^(rJ')) to 9'^^^. This fact, together the strong convergence 
of Pn in C{[0,T]; LP{n')) for p < 6™ — 3, implies that ^/l + + 6*°/^) converges to 
yrTp(l+r/2) in L^{n'x{0, T)). That means K^/'^{pn,en) strongly converges to 9) 
in L^{n' X (0,T)). Similarly, it follows that (1 4- p„)i/26'i''+'=+^^/^ converges strongly to 
in L^{Q'x{0, T)) due to the strong convergence of p„ and Therefore, 
Ko(/On,6'„)(l+/9„)i/2(l-h6l„)('^+"+i)/2 converges strongly to Ko{p, e){l + py/^{l + 9)'-''+''+^^/^ 
in L^{n' X (0, T)). On the other hand, we deduce from (g^) that (1 4- p„)^/^V(l + 
^'n)'"^^^^-'^^ is uniformly bounded in i^(f2 x (0, T)), hence weakly converges to some 
element uj in L'^{n x (0, T)). It also follows that V(l + 6'n)(°-'=+i)/2 jg uniformly bounded 
inL2(f7x (0,r)), and hence weakly converges to V(l 6')('^-=+i)/2 inL'^{n' x {0,T)). Due 
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to the strong convergence of p„ in L^{fl' x (0, T)), we deduce that lu = {1 + p)^/^V(l + 
^)(a-c+i)/2^ Finally, we write 

This, together with the strong convergence of Ko{pn,On){l + /0„)^/^(l + (?„)('^+^+^'/^ and 
the weak convergence of (1 + p„)^/^V(l + implies that k(p„, 6'„)V0„ converges 

to Ac(p, 6')V6' at least in the sense of distributions. 

For the viscous terms, y'^(p„)I?(u„) and A/|A(p„)|divu„ converges weakly to ^ p{p)D{u) 
and ■\/|A(p)|divu respectively in i^(fJ'x (0, T)), because of the hypothesis ()2.2p . the uniform 
bound on p„ in Lemma 13. 2[ and the uniform estimate (|4.ip . On the other hand, p,/ u„ 
strongly converges to p^^^u in L^{n' x (0,r)) and pn^^'^ \/ p{pn), as well as pn^^^ \/\{fin) 
converges strongly to p~^^^y/ p{p) and p~^/^ \{p) respectively in L°°([0, T]; L^(ri')). Hence 
^'„u„ converges to Vtu at least in the sense of distributions. 

Finally, we deal with the convergence of two nonlinear terms: (u„ x H„) x H„ and 
v{pn,On)^n X (V x H„). First, since H„ converges strongly to H in i''([0, T]; i^(ri')) for 
some p > 5 and some q > 5/2, u„ x H„ weakly converges to u x H in LP([0,T]; L'^{il')) 
for some p > 5/4 and q > 5/3 because u„ converges weakly to u in L^/^([0, T]; i^(rj')). 
From this, we can deduce that (u„ x H„) x H„ converges to (u x H) x H in the sense 
of distributions. Second, the strong convergence of p„, 0„, H„ and the hypothesis (|2.6|) 
imply that I'ipn, ^'ti)H„ converges strongly to i'{p, 9)11 in Lp{[0, T]; L'^{^1')) for some p > 5 
and some q > 5/2. By the weak convergence of H„ in L'^{[0,T]; H^{n)), one deduce that 
i^iPn, 6'ti)H„ X (V xH„) converges to i'{p, 6')Hx (Vx H) at least in the sense of distributions. 

Therefore, the energy conservation equation (|l.lcp holds at least in the sense of distri- 
butions. 

Step 4-' Convergence of the magnetic field equation. Similar to the argument in Step 
3, we can show that I'ipn, ^'n)V x H„ converges weakly to ^{p, 9)W x H in L^{Vl' x (0, T)). 
Also the strong convergence of H„ in i^([0, T]]L'^{^')) for some p > 5 and some q > 5/2 
and the weak convergence of u„ in L^/'^{[0,T];L^{fl')), imply that u„ x H„ converges to 
u x H at least in the sense of distributions. Hence, the induction equation holds at least 
in the sense of distributions. 

The proof is complete. 
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